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Zusammenfassung 


Es  werden  die  unteren  Eigenfrequenzen  einer  Kreisplatte  mit  gemischten  Rand- 
bedingungen  bestimmt.  Fur  eine  Kreisplatte,  deren  Rand  t eilweise  eingespannt , und 
gelenkig  gelagert  ist,  wurde  der  Verlauf  der  unteren  axialsymmetrischen  Eigenfre¬ 
quenzen  bestimmt.  Asymmetrische  Schwingungen  zeigen  jedoch  zwei  verschiedene 
Frequenzzweige  mit  jeweils  verschiedenen  angularen  Knotenlinien. 
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1  Introduction 


The  problem  of  bending  of  a  circular  plate  with  mixed  boundaries  has  been  treated  pre¬ 
viously  for  some  special  conditions.  A  more  general  problem  for  a  plate  partly  clamped 
and  simply  supported  was  formulated  in  [1],  where  forced  vibrations  due  to  harmonic  load 
perpendicular  to  the  plate  have  been  treated.  In  addition  the  plate  was  subjected  to  a 
compressive  load.  In  later  papers  [2-4]  a  variational  approach  was  applied  to  a  circular 
plate  partially  clamped  and  partially  simply  supported.  One  method  is  based  upon  two 
perturbations,  i.e.  one  when  the  plate  is  clamped  all  around,  the  other  when  the  plate  is 
simply  supported.  The  first  perturbation  yielded  upper  bounds  for  the  eigenvalue,  while 
the  latter  presented  the  lower  bounds.  There  are,  however,  still  some  discrepancies  of 
results,  all  of  which  have  been  restricted  to  the  mixed  boundary  conditions  of  a  partly 
clamped  and  partly  simply  supported  system,  for  which  only  the  fundamental  axisym- 
metric  mode  has  been  determined. 

With  the  advent  of  very  efficient  high  speed  computers,  exhibiting  the  advantage  of  solving 
a  large  number  of  algebraic  equations  at  a  relatively  short  time,  we  have  proposed  another 
method  for  the  determination  of  the  eigenvalues  of  various  mixed  boundary  conditions, 
where  not  only  the  fundamental  natural  frequency  has  been  determined,  but  also  higher 
axi-  and  asymmetric  modes  of  the  plate  were  investigated.  In  addition  we  have  determined 
the  nodal  lines  of  the  natural  frequencies  for  the  asymmetric  modes. 

It  should  be  mentioned  that  the  method  is  not  only  restricted  to  the  determination  of 
eigenfrequencies,  but  may  also  be  applied  to  radially  loaded  plates,  the  response  to  har¬ 
monically  perpendicular  loads  on  the  plates,  as  well  as  the  buckling  of  them.  In  addition, 
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the  method  may  also  be  used  for  other  boundary  combinations,  as  well  as  for  rectangular 
plates  with  mixed  boundary  conditions. 


2  Basic  equations 

The  problem  of  the  determination  of  the  approximate  lower  natural  frequencies  of  a 
circular  plate  exhibiting  partially  mixed  boundary  conditions  along  its  periphery  may 
be  performed  with  a  semi-analytical  method  as  shown  below.  This  method  satisfying 
the  boundary  conditions  for  a  finite  number  of  points  at  the  periphery  r  =  a  may  be 
applied  to  a  large  variety  of  boundary  conditions.  We  shall  treat  here  only  the  mixed 
boundary  conditions  of  clamped,  simply  supported,  while  those  of  free,  guided,  clamped, 
simply  supported  and  elastically  supported  edges  of  various  peripheral  edge  ranges  shall 
be  treated  at  a  later  point. 

The  basic  equations  require  the  solution  of  the  equation  of  the  circular  plate 

r\ 2 

DAAw  +  Qh-^-  =  0,  (1) 


where  the  operator  is 

A-—  1  9  1  92 

dr2  ^  r  dr  r 2  dp2 

Equation  (1)  has  to  be  solved  with  the  appropriate  mixed  boundary  conditions. 
The  bending  and  twisting  moments  are  given  by 


Mr  =  —D 


d2w 
dr 2 


+  v 


1  dw 
r  dr  ^ 


1  d2w\ 
r2  dp2) 


Mv  =  —D 


1  dw  d2w  1  d2w 
r  dr  ^  dr2  ^  r2  dip2 


(2) 

(3) 
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(a)  clamped:  w  =  0  and  ^  =  0  at  r  =  a 

(b)  simply  supported:  w  —  0  and  Mr  =  0  at  r  —  a 

(c)  free:  Mr  —  0  and  Vr  =  0  at  r  =  a 

(d)  guided:  ^  =  0  and  Vr  =  0  at  r  =  a 

(e)  elastically  supported:  Mr  —  =  0  and  Vr  +  kw  =  0  at  r  =  a, 

where  K  is  the  distributed  stiffness, i.e.  moment/unit  length,  opposing  with  spiral  springs 
the  edge  rotation,  and  where  k  is  the  distributed  stiffness,  i.e.  force/unit  length,  opposing 


the  translational  motion  in  direction  w. 


The  solution  of  equation  (1)  yields  with  w(r,  <p,  t)  —  W(r,  p)etojt  the  expression 


w{rM  =  E  {44  (a0  +  cmim  (a0| 


cos  mip 


sm  rrup 


where  A2  =  Loa2\Jgh/D.  The  problem  of  the  determination  of  the  approximate  lower 
natural  frequencies  for  mixed  peripheral  edge  conditions  depends  on  the  range  of  the 
various  boundary  conditions  at  hand.  It  may  be  noticed  that  certain  combinations  of 
mixed  boundary  conditions  show  for  the  total  range  one  and  the  same  boundary  condition, 
therefore  reducing  the  effort  of  the  numerical  procedure.  This  may  be  noticed  to  be 
obvious  for  a  combination  of  the  following  boundary  conditions. 


The  two  mixed  boundary  conditions  clamped  -  simply  supported  exhibit  at  r  =  a  the 
same  boundary  condition  w  —  0  at  r  —  a.  For  a  clamped  -  guided  mixed  boundary  case 
the  total  range  0  <  <p  <  2tt  shows  dw/dr  =  0  at  r  =  a.  For  a  system  of  simply  supported 
and  free  mixed  boundary  case  the  total  range  0  <  ip  <  2n  has  to  satisfy  the  vanishing 
bending  moment  Mr  —  0  at  r  =  a,  while  for  a  mixed  boundary  condition  free  -  guided 
the  Kelvin- Kirchhoff  edge  reaction  Vr  =  0  at  r  =  a,  i.e. 


Vr  =  -D 

Vr  =  -D 


S(A»)  + 


(1-y)  d2 
r  drdip 


-0, 


d3w  1  d2w  1  dw  (2  —  v)  d3w 
dr3  r  dr2  r 2  dr  r2  drdtp2 


or 

(3  -  1/)  a2wl  = 
r3  dip2 


(10) 


for  the  total  range  0  <  ip  <  2tc.  It  may  be  mentioned  that  for  a  ’’pure”  boundary 
condition,  i.e.  the  vanishing  range  a  =  0,  the  well-known  values  A2  of  that  particular 
boundary  condition  case  are  obtained,  while  for  the  ’’pure”  boundary  condition  of  the 
other  range,  i.e.  that  at  a  =  27T,  the  values  A2  of  that  particular  boundary  condition 
case  are  obtained.  For  mixed  boundary  condition  cases  of  varying  a^O  and  a  ^  2n  the 
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approximate  A2-values  are  obtained  from  the  following  procedure  and  require  for  different 
modes  ra  and  n  careful  selection  of  the  number  of  points,  at  which  the  remaining  boundary 
conditions  have  to  be  satisfied  to  yield  acceptable  final  results.  More  shall  be  remarked 
in  the  numerical  solution  below. 

3  Method  of  solution 

Let  us  treat  a  plate  clamped  along  part  of  the  boundary  and  simply  supported  along  the 
remainder,  requiring  only  two  different  boundary  conditions  for  the  numerical  procedure 
(Fig.  1).  If  the  plate  is  clamped  in  the  range  0  <  ip  <  a  and  simply  supported  in  the 
range  a  <  <  2x,  the  plate  exhibits  for  the  total  range  0  <  tp  <  2n  a  vanishing  deflection 

w  =  0  at  r  =  a  and  exhibits  therefore  the  solution 

W(r,ip)='%2  J m  (A-)  -  (a-)  {Amcosrrup  +  Bm  sin  imp}  .  (11) 

m=o  L  V  <*/.  MA)  V  aj 

This  satisfies  W(a,  y>)  —  0  in  the  total  range  0  <  ip  <  2n.  If  a  =  2i r,  i.e.  the  plate  is 
totally  clamped  the  second  boundary  condition  dw/dr  =  0  at  r  =  a  yields 

4(A)  -  =  0  (12) 

and  we  obtain  for  the  determination  of  Xmn  the  vanishing  determinant 

Jm(A)  Im(A) 

=  0  (13) 

J„(A)  t(A) 

which  yields  the  squares  of  the  eigenvalues  A^2  for  a  clamped  plate  given  in  table  1. 

If  the  plate  is  totally  simply  supported,  the  value  a  =  0  and  equation  (11)  has  to  satisfy 

in  addition  the  vanishing  bending  moment  Mr  =  0  at  the  edge  r  —  a.  This  yields  with 
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(2) 


Jm(A)  +  j4(A)]  Im(A)  -  [t(A)  +  jtw]  =  0  (14) 

and  the  determination  of  A  requires  the  solution  of  the  determinant 

Jm(A)  WA) 

=  0,  . (15) 

1»+M  CW  +  &.W 

which  yields  for  v  =  0.3  the  squares  of  the  eigenvalues  A^2  for  a  simply  supported  plate 
given  in  table  2. 


For  a  plate  of  partly  clamped  boundary  in  the  range  0  <  <p  <  a  and  a  simply  supported 
boundary  in  the  range  a  <  tp  <2ir  (Fig.  1)  we  have  to  satisfy  dw/dr  —  0  at  r  =  a,  i.e. 

°°  r  j  m 

E  Jm  (A)  ~  (A)  (Am  cos  rmp  +  Bm  sin  rmp )  =  0  in  the  range  0  <  (p  <  a 

m=0  L  bn  (A)  J 

(16) 

and  Mr  =  0,  i.e. 

E  {  Jm  (A)  +  (A)  -  C  (A)  +  jC  (A)  |  (Am  cos  rrup  +  Bm  sin  rmp)  =  0 

in  the  range  a  <  <p  <  2rr.  (17) 


The  equations  (16)  and  (17)  have  to  be  satisfied  at  a  chosen  number  of  points  in  each 
range.  If  <p  —  ani/(Ni  +  1)  with  n\  =  1,  2,  . . . ,  N\  in  the  range  0  <  <p  <  a  and 
(p  =  a  +  (2rr  —  a)n2/iV2  with  n2  =  0,  1,  . . . ,  iV2  in  the  range  a  <  <p  <  2n,  the  equations 
(16)  and  (17)  read  then  {N\  +  i V2  even) 


(Ni+m/ 2  r  T  Q'i 

E  J'mW-fm1 

m=0  bnfAJ 


»(A)  { 


Am  COS 


ma  n\ 


iVi  + 


n  .  n  ■  mani  \  i 

lJ+B”*sln(^TT)} 


for  n\ 


8 


and 


(N1+N2)/2  f 


ro=0  (  L 


E 

x  ^Am  cos 


v 


J"(A)  +  ^4(A)j 

(27r  -  a)ri2 


Jm(A) 


V , 


a  + 


IVo 


I»(A)  L 

+  Bm  sin  fra 


C(A)  +  X4(A) 


x 


a  + 


(27r  -  a)ri2 


iVo 


for  n2  =  0,  1,  . . . ,  jV2. 


(19) 


The  equation  (18)  represents  N\  homogeneous  algebraic  equations  in  the  unknowns 
4o,  Au  . . . ,  A{Ni+N2)/ 2  and  i?i,  B2>  B{Nl+N2)/ 2,  while  equation  (19)  yields  1V2  +  1 

equations  in  those  constants.  The  vanishing  coefficient  determinant  represents  the  equa¬ 
tion  for  the  determination  of  the  lower  approximate  eigenvalues  A.  It  yields  the  lower 
values  \mn  for  a  given  magnitude  of  a  in  ascending  order  of  which  for  a  given  mode  m 
the  values  A2  shall  be  located  between  those  of  the  totally  clamped  case  and  those  of  the 
totally  simply  supported  case. 


4  Numerical  evaluations  and  conclusions 

The  above  procedure  has  been  numerically  evaluated  for  the  lower  modes  of  a  circular 
plate,  i.e.  for  m  =  0,  1  and  2  and  n  =  1,  2.  Figure  2  represents  A2at  =  ua2^Jgh/D 
as  a  function  of  the  angle  a/n.  For  a  =  0  we  exhibit  the  plate  with  a  purely  simply 
supported  boundary,  while  for  a  =  2n  the  boundary  is  in  a  purely  clamped  state.  These 
values  are  indicated  in  the  graph  as  o  for  the  simply  supported  boundary  and  as  0  for 
the  clamped  boundary.  The  influence  of  the  Poisson  ratio  u,  which  appears  only  in  the 
simply  supported  boundary  is  indicated  by  the  *-sign  where  the  upper  star  represents 
the  A-value  for  v  =  0.5  and  the  lower  one  that  of  v  —  0.2.  The  square  of  the  eigenvalue 
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^01  increases  from  4.9351  to  10.2158  and  exhibits  a  slight  flexion  close  to  a  =  it.  The 
here  presented  numerical  results  were  compared  with  those  of  reference  [3]  and  [4],  where 
only  the  mode  m  =  0,  n  =  1  has  been  treated  and  presented.  The  value  presented  in 
Fig.  2  for  m  =  0  and  n  =  1  show  a  very  close  proximity  to  that  given  in  [3].  The 
numerical  evaluation  in  our  treatment  required  for  Ni  +  N2  =  150,  in  which  N\  and  N2 
vary  according  to  the  magnitude  of  a.  The  magnitudes  of  N\  and  i\T2  are  chosen  such  that 
an  increase  of  these  magnitudes  do  no  longer  affect  the  accuracy  of  the  plotted  A2-values. 
If  a  is  small  the  portion  of  the  boundary  being  clamped  needs  only  a  small  value  for  N\ , 
while  that  of  the  boundary  being  simply  supported  requires  a  large  number  N2  for  the 
numerical  procedure.  As  the  clamped  portion  a  of  the  boundary  increases,  so  does  the 

t 

number  N\,  while  N2  decreases,  such  that  the  sum  N\  +  N2  =  150,  as  indicated  in  Fig. 
2.  Treating  the  axisymmetric  mode  m  =  0  further  yields  the  results  shown  as  the  curve 
m  =  0,  n  =  2  (second  axisymmetric  mode),  starting  for  a  totally  simply  supported  plate 
at  Aq2Ss)  =  29.72  and  reaches  for  totally  clamped  plate  the  value  A^  =  39.77  (see  also 
Table  1  and  2).  In  the  case  of  the  eigenvalue  A22  we  notice  increased  flexion  as  a  increases. 
The  investigation  of  the  asymmetric  mode  m  =  1  and  m  =  2  exhibits  for  the  first  modal 
number  n  =  1  two  values  A12  and  A2i  in  the  range  0  <  a  <  27r.  For  m  =  1  and  n  —  1 
the  value  of  A2X  yields  for  a  totally  supported  plate  the  magnitude  A^ss')  =  13.8982  and 
for  a  completely  clamped  plate  the  value  A^  =  21.26.  In  the  mixed  boundary  range 
a  ^  0,  2n  we  have  to  distinguish  two  modes,  which  have  been  indicated  by  Q)  and  ©. 
The  mode  Q)  exhibits  in  the  range  0  <  a  <  7r  eigenvalues  Af 2  being  smaller  than  those 
of  the  branch  while  above  a  >  n,  i.e.  the  case,  where  more  than  half  of  the  plate  is 
clamped  at  its  boundary,  the  eigenvalue  X2n  is  larger  than  that  of  the  branch  Those 
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two  branches  caused  by  the  mixing  of  boundaries  shall  not  only  exhibit  different  natural 
frequencies,  but  also  different  nodal  lines,  as  we  shall  see  later.  Treating  now  the  mode 
m  =  2,  n  =  1,  which  shall  start  for  a  completely  simply  supported  plate  with  the  square  of 
the  eigenvalue,  i.e.  A^f*^  =  25.61  and  reaches  for  the  totally  clamped  plate  the  magnitude 
=  34.877,  exhibits,  also  two  branches  (I) and  (2) with  increased  fluctuations.  In  Fig.  3 
we  present  the  fundamental  axisymmetric  mode  to  =  0,  n  =  1  again,  where  N1  +  N2  —  160 
and  v  was  chosen  to  be  v  —  0.2,  0.3  and  0.5.  If  the  Poisson  ratio  v  is  larger  the  natural 
frequency  is  larger,  the  magnitude  of  which  decreases  as  a  increases.  A  similar  behavior 
is  found  for  a  Poisson  ratio  v  being  smaller.  Figure  4  shows  the  numerical  results  for  the 
mode  to  =  1,  n  —  2,  which  exhibits  again  two  branches  for  the  natural  frequency  of  the 
plate  for  a  ^  0  and  2n.  Both  curves  start  for  a  =  0  at  =  48.48  and  end  for  a  —  2ir 
at  A5c)  =  60.83. 

In  Fig.  5  we  represent  the  nodal  lines  for  the  branches  Q  and  Q  of  the  mode  to  = 
n  =  1  with  the  angle  a  as  parameter  in  degrees.  Behind  each  a  the  magnitude  of  the 
dimensionless  natural  frequency  A2  =  coa2yJgh/D  =  ua?  s]V2q{1  -  v)/EIh 2  is  indicated. 
The  drawing  of  the  node  lines  was  done  as  follows.  For  chosen  a-value,  mode  numbers 
to,  n  and  branch  the  A2 -value,  i.e.  the  eigenvalue,  and  the  coefficients  Am,  Bm,  i.e.  the 
eigenvector,  are  computed  from  (18)  and  (19).  Introducing  these  results  into  (11)  and 
demanding  W(r,<p)  —  0  yields  for  chosen  r/a  with  (0  <  rja  <  1)  the  appropriate  angle 
ip  =  </?(r).  With  increasing  angle  a,  i.e.  increasing  portion  of  the  clamped  boundary 
the  nodal  line  turns  in  counter-clockwise  direction,  as  indicated.  For  a  totally  simply 
supported  plate  the  nodal  line  is  located  at  ip  =  0.  This  is  also  valid  for  a  completely 
clamped  plate,  while  for  a  plate  of  which  half  of  the  boundary  is  clamped  and  the  other 
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half  simply  supported,  the  nodal  line  is  located  at  <p  —  7t/2,  i.e.  perpendicular  to  that  of 
the  pure  boundary  cases.  If  only  one  quarter  of  the  boundary  is  clamped,  i.e.  a  =  7r/2, 
then  the  nodal  line  is  at  <p  =  7r/4.  The  nodal  lines  for  branch  (2)  of  the  mode  ro  =  n-  l 
are  also  presented  in  Fig.  5.  These  nodal  lines  appear  as  soon  as  a  ^  0,  2 it.  For  a  small 
clamped  part  of  the  boundary  of  a  —  10°,  for  which  =  14.80  the  nodal  iine  exhibits  a 
small  curvature  and  does  no  longer  pass  through  the  center  of  the  plate.  As  a  increases  the 
curvature  of  the  nodal  line  increases  and  rotates  as  a  whole  in  counter-clockwise  direction, 
exhibiting  increased  curvature  especially  at  the  area  close  to  the  boundary.  While  for  the 
branch  0  and  a  =  it  the  nodal  line  is  a  straight  line  4>  =  7t/2,  that  of  the  branch  (2)  shall 
be  located  as  a  curved  line  in  the  lower  half  (i.e.  the  simply  supported  part)  of  the  plate 
and  shall  be  below  the  center  of  the  plate  with  a  larger  inward  curvature  close  to  the 
boundary.  This  curvature  will  increase  strongly  with  the  increase  of  the  clamped  part  a. 

For  the  branches  Q  and  ©  of  the  mode  m  =  2,  n  —  1  the  nodal  lines  are  presented  in 
Fig.  6.  If  the  clamped  portion  of  the  boundary  is  a  small  a  =  10°  the  nodal  lines  of 
branch  ©are  the  two  straight  and  perpendicular  lines  ©as  indicated  in  Fig.  6a,  rotated 
in  a  counter-clockwise  direction  by  5°.  The  nodal  lines  associated  with  branch  ©are  the 
exhibited  hyperbola-like  curves.  A  further  increase  of  a,  meaning  a  larger  portion  of  the 
clamped  boundary  yields  a  counter-clockwise  rotation  of  the  nodal  curves  ©  and  ©  It 
may  be  noticed  that  perpendicular  nodal  line  of  the  branch  ©shifts  from  the  center  of  the 
plate  and  exhibits  slight  curvature  at  the  boundary  of  the  plate  (see  Fig.  6c  for  a  =  120°). 
In  addition  the  nodal  lines  of  branch  ©show  at  those  nodal  curves,  close  to  the  clamped 
boundary  increasing  downward  curvatures,  as  a  increases.  For  larger  a- values  the  nodal 
lines  are  presented  in  Figs.  6d  through  6g. 
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eigenvalue 


13 


References 


[1]  Nowacki,  W.;  Olesiak,  Z.:  Vibration,  buckling,  and  bending  of  a  circular  plate 
clamped  along  part  of  its  periphery  and  simply  supported  on  the  remaining  part. 

.  Bull.  Acad.  Pol.  Sci.,  cl.  IV,  4(4)  (1956)  247-258 

[2]  Bartlett,  C.;  Noble,  B.A.:  A  variational  method  for  the  solution  of  eigenvalue  prob¬ 
lems  involving  mixed  boundary  condition.  Appl.  Sci.  Res.  B  9  (1962)  403-419 

[3]  Bartlett,  C.:  The  vibration  and  buckling  of  a  circular  plate  clamped  on  part  of  its 
boundary  and' simply  supported  on  the  remainder.  Quart.  Joum.  Mech.  and  Applied 
Math.  16  (1963)  431-440 

[4]  Noble,  B.:  The  vibration  and  buckling  of  a  circular  plate  clamped  on  part  of  its 
boundary  and  simply  supported  on  the  remainder.  Proc.  9th  Midwest.  Conf.  Solid 
and  Fluid  Mech.  (1965)  141-146 


14 


List  of  Tables 


Table  1:  Squares  of  the  eigenvalues  for  a  clamped  plate 

Table  2:  Squares  of  the  eigenvalues  A^2  for  a  simply  supported  plate  for  Poisson  ratio 

v  =  0.3  ■'  . . ' . .  '  ' 


15 


n  \  m 

0 

1 

2 

3 

1 

10.2158 

21.2604 

34.8770 

51.0300 

2 

39.7711 

60.8287 

84.5826 

111.0214 

m 

89.1041 

153.8151 

190.3038 

4 

158.1842 

199.0534 

242.7206 

289.1799 

Table  1:  A^2  for  a  clamped 


n  \  m 

0 

1 

2 

3 

1 

4.9351 

13.8982 

25.6133 

39.9573 

2 

29.7200 

48.4789 

70.1170 

94.5490 

3 

74.1561 

102.7733 

134.2978 

168.6749 

4 

138.3181 

176.8012 

218.2026 

262.4847 

Table  2:  A^2  for  a  simply  supported  plate  (u  =  0.3) 
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